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Abstract 

The number of hydrogen-like atoms produced when heavy nuclei collide is esti- 
mated for central collisions at the Relativistic Heavy Ion Collider using the sudden 
approximation of Baym et al. As first suggested by Schwartz, a simultaneous mea- 
surement of the hydrogen and hadron spectra will allow an inference of the electron 
or muon spectra at low momentum where a direct experimental measurement is not 
feasible. 
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The production rate of lepton pairs is a rapidly increasing function of temperature 
and so has long been considered a good probe of the initial high energy density phase of 
ultrarelativistic nuclear collisions The experimental detection of such direct leptons 
is a problem in the sub-GeV range of transverse momentum due to the large number of 
charged hadrons produced and the need to disentangle direct leptons from those arising 
from hadron decays. But this is just the kinematic range characterizing a quark-gluon 
plasma at a temperature of 200 to 500 MeV. 

Schwartz || proposed to measure the distribution of atoms formed by the binding of 
a directly produced lepton to one of the charged hadrons emerging from the final state of 
the nuclear collision. A measurement of the charged hadrons and of the atoms, together 
with a theoretical calculation relating the distributions of the three particle species, would 
then imply the spectrum of leptons. The beauty of the idea lies in the fact that nearly 
all indirectly produced leptons arise from the decay of hadrons, and these decays occur 
too long after the collision to allow an atom to be formed. Of course, one still cannot tell 
whether the leptons were produced in quark-gluon plasma or in hadronic matter, but this 
is another issue. 

Five years ago Baym, Friedman, Hughes and Jacak calculated the relationship among 
the spectra of the atoms and of the charged hadrons and leptons which comprise them 
||. The formula reads: 

dyd 2 p ±tatom md dyd 2 p ±:h dyd 2 p ±j i ' 

Here £(3) = 1.202... and m re d is the reduced mass of the hadron and lepton making up the 
atom. Since the binding energy is so small it is an excellent approximation to evaluate 
the hadron and lepton rapidities at the same rapidity as that of the atom, and to equate 
their transverse velocities as well: p_i_, a tom/ m atom = P±,h/ m h = P_L,i/ m i- This formula 
is based on the sudden approximation which simply asks for the overlap of the outgoing 
wave functions of the hadron and the lepton with their hydrogenic state. The sudden 
approximation is valid because these particles are formed in a nuclear volume which is 
extremely small in comparison to the size of the hydrogen atom and over a time interval 
which is extremely small in comparison to the Bohr period. The specific focus of Baym 
et al. was on 7r-/i atoms. Here we shall be interested in p-e, p-/z, 7r-e and 7r-/i atoms. 
Our essential contribution is to estimate dN/dyd 2 p± for the leptons, protons and pions 
in the relevant range of transverse momentum, and from these to estimate the number of 
hydrogenic atoms to be formed in central Au+Au collisions at the Relativistic Heavy Ion 
Collider (RHIC). 

First we estimate the number of leptons produced in the quark-gluon plasma phase. 
The reaction rate for the process q + q — > l + + l~ is: 

p 19 [ d3pi d3p2 d3p+ ^Ez 

q J 2E l {2nf 2£ 2 (27r) 3 2£+(27r) 3 2£_(27r) 3 

x f FB (Ei)fFB(E 2 )\M\ 2 (2n) i 6\ Pl +p 2 -p--p + ), (2) 
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where the 12 arises from three colors and four possible spin states of the colliding quarks, 
/fd is the Fermi-Dirac distribution, and Ai is the matrix element for the reaction. Ap- 
proximating /fd(-E) ~ e~ E l T the momentum distribution for negatively charged leptons 
becomes [ffl: 



d 3 R ^ e 2 e 2 T 



(Pp_ ^ (2vr 

' q=u,a,s v 



V 6 6g e~ g - /T /" ds In (l + e~ s/4B - T 



x [l + 2(m 2 q + m 2 )/s + 4m 2 q m 2 /s 2 ] ^(1 - Am 2 /s)(l - Am 2 J s) , (3) 

where s = (p± + P2) 2 - The expression for positively charged leptons is the same. If the 
masses of the leptons and the quarks can be neglected this simplifies to 

For the range of transverse momentum of interest in this context the masses of the muon 
and strange quark cannot be neglected. 

To obtain the total number emitted we integrate over the space-time volume according 
to Bjorken's model ||. For central collisions: 

dN ™ fe q \ 2 3a 2 , ^ A 2 dr] [ T ° dT 



dyd 2 pT 



q=u,d,s 



x [l + 2(mJ + m 2 )/s + 4mjmf / s 2 ^ (1 - Am 2 /s)(l - 4m 2 q /s) . (5) 

Here R? is the nuclear radius, y is the momentum space rapidity, r\ is the position space 
rapidity, T is the temperature when the plasma is first considered to be thermalized, and 
T c is the critical or phase transition temperature. In the Bjorken model the temperature 
drops with proper time r according to 

T(r) = j T . (6) 



Finally E = m± cosh(?/ — 77) and mj_ = ymf + pj_ where the z-axis is the beam axis. 

In general the integrals must be done numerically. However, if the masses can be 
neglected then the above simplifies to: 

dN 2a 2 R 2 T (r T 3 ) 2 

dyd 2 pT vr 3 pj, 

-(p T /T c fKi(p T /T c ) - 2(p T /T c ) 2 K 2 ( PT /T c )} . (7) 



(pT/TofK^pr/To) + 2(p T /T ) 2 K 2 (p T /T ) 



If one is interested in hydrogen-like atoms with a transverse momentum of a few GeV/c 
then the transverse momentum of the lepton must have been 



mi 

P±,l = Rl_,atom 

^atom 



which is just a few MeV/c for electrons and a few hundred MeV for muons. For electrons 
the above formula simplifies even more since their transverse momentum is always much 
less than the temperature. 



dN a Rjn 2 



dyd 2 pT 27r 3 



T \ 4 ln A.20T C \ _ ln /l.207y 



T C J \ P± J \ P± 



(9) 



This invariant distribution diverges logarithmically at small transverse momentum. 

Our estimates of the lepton distributions also include those coming from the mixed 
phase as the system goes through a first-order phase transition. The contribution from 
the quark-gluon plasma phase at the phase transition temperature T c is given by the 
distribution (|3|) multiplied by the volume fraction occupied by the plasma, / p i a sma = 
(rr c /r — 1) /(r — 1) plus the corresponding rate in the hadronic phase multiplied by the 
hadronic volume fraction /had = 1 — /plasma- Here r is the ratio of the number of degrees 
of freedom in the two phases. The hadronic rate is obtained from the annihilation process 
7r + + 7r~ — ► / + + l~. These calculations are standard and well-known [Q. 

The canonical picture of central collisions of gold nuclei at RHIC is that the central 
rapidity region will be almost baryon free ||. However, central collisions of lead nuclei at 
the SPS |I| and extrapolations of nucleon collisions such as LEXUS suggest that the 
baryon rapidity distribution may be roughly flat. Furthermore, initial temperatures at 
RHIC may be as high as 500 MeV [0] and this will be reflected in the final transverse mass 
distribution of the outgoing protons since high temperatures eventually get converted to 
transverse flow. Indeed, inverse slopes for protons for central lead collisions at the SPS 
already reach 300 MeV 0. Therefore a not unreasonable estimate for protons at RHIC 
is to assume a flat rapidity distribution times an exponential falloff in transverse mass. 

Z exp [{pip - m ±p )/T p ] (10) 



dyd 2 p± 2iry (m p + T p )T p 

Z is the charge of a single nucleus, 2y is the rapidity gap between projectile and target 
nuclei, m± is the proton's transverse mass, and T p is the proton effective temperature 
(not really a temperature, just an inverse slope). Similarly, we assume the charged-pion 
distribution to be: 

Z exp [(m^ - m± n )/T^\ , (11) 



dyd 2 p± dy 2ir(m w + T 7T )T 7T 

with the pion rapidity distribution dN n /dy = 1000. 

To get a rough idea of the numbers assume: T = 3T C = T p = T n = 480 MeV, Z = 79, 
Rt = 7 fm for a central collision at the maximum RHIC energy of 100 GeV per nucleon 
per beam, r = 19/3 for three quark flavors and a hadron gas with = 7.5 effective 
degrees of freedom 0. 

In figures 1 and 2 we plot the transverse momentum distributions of p-e, 7r-e, p-fi and 
7r-/x atoms for a central Au+Au collision at RHIC. The trends are easily understood on 
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the basis of eq. (1). The muonic atoms dominate the electronic ones because they have 
a greater reduced mass, and the pionic ones dominate the protonic ones because charged 
pions are much more abundant than protons. 

In figures 3 and 4 we plot the number of hydrogen-like atoms expected per unit of 
rapidity per day for transverse momentum bigger than p± >m i n for Au+Au collisions with 
an impact parameter of 1 fm or less. A beam luminosity of 2 x 10 26 /cm 2 sec is assumed. 
For the most abundant species, ir-/i atoms, we estimate about 1000 per unit of rapidity 
per day for transverse momenta larger than 1 GeV/c. 

In conclusion, we have reinvestigated the rates for the production of hydrogen-like 
atoms at RHIC. The results are quite promising for their experimental detection. It 
remains to be seen whether an efficient detector can be designed to observe them. 
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Figure Captions 



Figure 1: The transverse momentum distributions for electronic atoms for a central 
Au+Au collision at RHIC. 

Figure 2: The transverse momentum distributions for muonic atoms for a central Au+Au 
collision at RHIC. 

Figure 3: The number of electronic atoms produced with a transverse momentum greater 
than the indicated value per unit of rapidity per day at RHIC. These assume design lu- 
minosity and impact parameters less than 1 fm for Au+Au collisions. 

Figure 4: The number of muonic atoms produced with a transverse momentum greater 
than the indicated value per unit of rapidity per day at RHIC. These assume design lu- 
minosity and impact parameters less than 1 fm for Au+Au collisions. 



6 



CM 



CM 



-a 



10 



10 



> 10 

o 

' 10 

H 



-9 



-10 



-11 



10 



10 



10 



10 



-12 



-13 



-14 



-15 



10 



-16 




1 2 3 4 5 6 

P T , mi n [GeV/c] 



